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Abstract 


The particle density, and the magnetic field intensity and direc- 
tion, are calculated in corotating streams of the solar wind. It is 
assumed that the solar wind velocity is constant and radial and that its 
azimuthal variations are not too rapid. These assumptions are approxi- 
mately valid between 0.1 and 1 AU for many streams. In the absence of 
streaming, n decreases as r“^ and B/n is relatively constant, although 
modified by corotation. Kinematic factors change the variation of n 
somewhat. Their effect on B is similar, but also depends on the initial 
orientation of When our assumptions are valid, dynamic effects, 
which are considered briefly, will also change n and ^ but to a lesser 

-2 

degree. By averaging over a typical stream, it is found that ~ r > 

whereas '^B 0 > does not vary in a simple way. Changes of field direction 
may be very large, depending on the initial angle; but when the initial 
angle at 0.1 AU is such that the base of the field line corotates with 
the sun the spiral angle is the preferred direction at 1 AU, 
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1. Introduction 


The gas-dynamic properties of interplanetary streams have been 
extensively studied and were reviewed by Hundhausen (1972) and Burlaga 
(1974) , but relatively little has been written about the theory of 
magnetic fields in streams. Sakurai (I971) considered a model of 
stationary, co-rotating streams which includes magnetic fields, B, but 
neglects the reaction of B on the velocity V. Matsuda and Sakurai 
(1972) and Urch (I972) extended this work to include, to first approxi- 
mation, the modifications of V induced by the pressure gradients that 
are produced by the steepening of a speed profile. Exact numerical 
solutions for two-dimensional, stationary corotating streams were 
obtained by Nakagawa and Wellck (1973) who introduced a temperature 
variation as well as a speed variation. All of these models consider 
an inner boundary at 30 and neglect the primary acceleration 

mechanism. We shall follow a similar approach; we consider that the 
streams begin as speed variations at 0.1 AU ~ 20 and determine the 
effects of the steepening of the stream on the magnetic field and 
density between 20 and 1 AU. 

The results of the stream-magnetism models mentioned above are all 
basically the same. They predict an enhancement of B in the leading 
part of the stream, a rarefaction in the trailing part and some per- 
turbations of the magnetic field direction from the spiral angle. 

However, all of these models assiane that B is nearly radial at 30 I^, 
and thus fail to consider the large fluctuations of field direction about 
the spiral angle which are a general feature near 1 AU. Schatten 
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(1972) has briefly discussed the effect of speed gradients on the 
direction of B and showed that the effect can be appreciable, but he 
considered only a few isolated events rather than stream profiles, and 
he did not examine magnetic field intensity variations. 

In this paper we systematically examine the effects of the radial 
velocity profile in corotating streams on the magnetic fields between 

J TT 

20 and 1 AU, using the kinematic approximation, ^ =0 (i.e,, V = 
constant for any given volume element) and a variety of B configurations 
on the inner boundary. Our approximation is valid for typical streams 
of moderate amplitude . It does not exactly 

describe very steep streams, but it allows one to examine a greater 
variety of conditions than one could explore with the relatively few 
solutions that one can compute in practice for more complicated dynamical 
models. Thus, the kinematic approximation provides considerable insight, 
which is our main objective. In any case, the zeroth order effects of 
the streams on B are kinematic effects. For example, the magnetic field 
compression obtained in the dynamical models referenced above, is a 
kinematic effect; dynamical changes in V are a consequence of this com- 
pression rather than the cause of it, 

2. Basic Equations 
a) Approach 

In practice, it is not possible to measure functions such as B (x, 
y, z, t) and V (x, y, z, t) at every point in the solar wind. Rather, 
one measures functions of time, at one or two points in space, which to 
first approximation are the result of the passage of a continuous series 
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of volume elements moving radially past the observer. Thus, it is 
reasonable to adopt a Langrangian point of view in which one follows 
volume elements moving on a radial line which joins the observer and 
the sun. Knowing the properties of each volume element as it moves from 
the sun to the observer, one can construct the time profiles that the 
observer sees as a result of the passage of a continuous train of such 
volume elements. This approach is used below, 

b) Velocity Gradients 

As will be shown later, the changes in the properties of a volume 

element as it moves radially from (r) to (r + Ar) depend on the local 

AV 

gradients ^ and 

reference point on the sun (See Figure 1). We obtain expressions for 
these quantities as follows. 

Let us assvime that at a distance rQ and an angle 8 in the frame 
rotating with the sun there is a stationary solar wind speed distribu- 
tion Vq( 8) a number density distribution n^(8) and a magnetic field 
distribution ^(8) (See Figure 1). As these distributions rotate past 
the observer-sun line, they generate time profiles which fom a stream 
on that line. Let us define a(8) in the corotating frame as ^ = a(8). 
Now let us consider a fixed frame of reference, in which the sun rotates 
with angular velocity Q =flz. We assume that the solar wind velocity 
in the fixed frame is everywhere radial, as a first approximation. 

Plasma emitted at time tj_ will reach r^^ at time Tj = t-j_ + 

At a time t 2 , the sun will have rotated by an amount | A8l = fl (t^-t^), 

A volume element emitted then will reach at a time 


^ at r, where 0 is the azimuthal angle from a 
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T 2 = t 2 + (^2"^o^/^ 2* sufficiently 

small, V 2 -V^ = a|5|=afl (t 2 "t^). Now consider volume elements arriving 
at rj^ and at the same time T; for these, t^_ + (rjL-r 2 )/Vj_ = tg + {^ 2 - 
t^)/M 2 - Setting r^ = r^ + Ar, a little algebra shows that 

. 0, JT-ViVt. 

VtCo.>n.(v»»V'^'43 


which becomes, in the limit Ar 


where 


Ic ^ 


with 


qA g-lp 
v'/ Ro 




Note that fi >0 and a>0 in the "rise" of a stream, so /i >1 in a region 
where V increases with time. One can similarly calciilate — as follows. 
Consider two volume elements moving with speeds V 2 and Vj_ at r at = T^, 
The difference between the departure times from r^ is tg-t^ = (r-r^) 

- V"^) , during which time the sun rotates through an angle A0^ = H 
(tg-t^). Considering the geometry at r^, one finds that + a(0 + 

A0I;.) . Collecting terras and taking the limits A0 *>0, At ►O, one 
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obtains 



Note that equations (1), (3), and (4) are valid for all latitudes, 9^. 


c) Density Variations 

We are now in a position to calculate the density variations in a 
volume element along its trajectory. Starting with the equation of 
continuity 


f tt' 


vrhere 9 = nm is the mass density, and assuming that V is constant along 
a trajectory, so that ^ ^ obtain 

^ ^ ^ y dvi . ( 6 ) 

f '>1 ^ ■ 0V t. 

Substituting (2) into this expression and integrating from r^^ to r 
gives 



(7) 


Our first result, then, is that the density does not decrease with 
distance exactly as r“^, unless a = 0 (i.e. no stream); a kinematic 
correction factor is necessary, equal to ft which is given by (3). The 
conditions under which (7) is valid are discussed in the Appendix. One 
can, however, state immediately that e:^n if the density at the source 
is constant, independent of S, a volume element in which the speed is 
increasing (a > o) will have a higher density at r > than an element 
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in which a = 0, since fi> 1 for a > 0. The variation of ft with a/v, for 
observers at 0,3 AU, 1,0 AU, and 1.5 AU, is shown in Figure S. The 
variation of with r is also shown in Figure 2, for a few values of a/v. 
The compression or rarefaction of n implied by the variation of fi has 
been noted by others in gas-dynamical models of streams (for references, 
see Hundhausen, 19T2), but it is interesting to note that this is a 
kinematic effect. The dynamical effects tend to reduce the amount of 
compression or rarefaction; dynamical effects are the result of 
the density change rather than its cause. 


d) Magnetic Field Variations 

We wish to find the variations of B(r) in a volume element along 
its trajectory. We start with the equation 




(8) 


which follows from one of Maxwell's equations with the frozen-field 
condition E = -V x B. Using a well-known identity for V x (V x B) 
with V • B = 0 and the definition of the convective derivative, d/dt = 


+ X * Sj we obtain 


^ =(?-v)y-6(v.y) 

... .-„-.„nt to consider the 
Using (5) and ( 9 ) one finds that 

i (i\ = J-fB VW 


(9) 


It is convenient to consider the Arariations of B/n rather than B directly. 

( 10 ) 


( 11 ) 
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This is known as Wal^n's equation (WalSn, 1946). 

The LHS of (ll) gives the change in B/n in a volume element which 
moves radially at constant speed relative to a fixed frame, as seen by 
an observer in the fixed frame, let us consider a spherically symmetric 

A 1 ^ • 

coordinate system in which r points radially , away from the sun, p is. in 
the direction of motion of the planets, and ^ is directed along ^ x r 
(See Figure 1, or see Stratton, 1941, p. 52). The co-latitude is 9 
(sin 0 = 1 in the ecliptic plane). 

A general solution of Walin’ s equation is found in many textbooks 


on plasma dynamics, viz.. 


l.j. (6. V)/?, 

-n. \ ~ - ' ~ 


where R is the displacement vector of a volume element (We use g instead 
of r in this discussion to avoid confusion between Lagrangian and Eulerian 


coordinates). However, this is often given in an incorrect form or else 
the meaning is not fully explained. One of the clearest explanations of 
this equation is found in Batchelor (1970). Eq, (12) is the form of the 
solution given by Boyd and Sanderson (I969). Let the 0, 0 components of 
B/n be denoted as 

respectively and let f^ and gg be the values at the inner boundary. 
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other derivatives are obtained by differentiating the equation R = + 

•o 

obtains 


V (ro, 0o) t at a given t, noting that^^ = - J 


One 



(15) 


^ V y- 

^ - -' -alR.?.) (,e) 

Inserting these expressions into (l4) replacing R by r, using (13) gives 


i.e 


• 5 


' \ yit, u, V* »»*, 'J»6vhQ/ 



(17) 


(18) 




(19) 


AtH & 

Since there is confusion about the meaning and derivation of (12), 
it is instructive to derive the solution of (ll) written in the particu- 
lar form appropriate for our problem. We can set dt = dr/v, because we 
are following a volume element, whose speed is constant in our picture. 
Assuming that Bq=0 and that V = Vr, (11) gives 

-hj -u 'JV- 

U:r A 


( 20 ) 


( 21 ) 
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Using (2) and (i<-) , Equation (21) can be integrated immediately, giving 
(l8). Substituting this result into (20) and integrating gives (19) • 


One can write (17) in the form 



When a = 0, g- = 0 and u = 1, and the components of B are just B 
(r^/r)^ and (r^/r). If, furthermore, B^ = B^ /R^, one has the 

result of Parker (1958) for a field line whose base corotates with the 
sun, assuming symmetric flow. In the case of a corotating stream, where 
a 0, B0 is enhanced by the factor pt, and B^ is augmented by an amovuit 
proportional to a and to B0 . When the field is initially radial, B = B^^ 
(r^j/r) ; in our approximation, the streaming does not modify B in this case. 
For the rather unlikely case of an initially azimuthal field, a radial 
component whc^e magnitude can be calculated from (22) is produced by 
the velocity shear in the stream. 
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3. STREAMS AT 1 AU 


Let us consider an observer at a fixed point at 1 AU and ask what 
he should measure as a corotating stream moves past him. For simplicity, 
let us begin by considering a stream whose speed profile at 0.1 AU is 
V Vq + AV cos 48, where = 400 km/s and AV = 75 km/(sec<,rad. ) Such 
a profile corresponds to a minimum speed of 325 km/sec and a maximum 
speed of 475 toti/sec, and to 4 streams per solar rotation. This approxi- 
mates the characteristics of the streams observed by IMP 1 (Wilcox and 
Ness, 1965) S’lid a stream which Hundhausen (1972) chose to fit with his 
gas-dynamic model. It may be regarded as a representative stream (e.g. 
see Hundhausen, 1972; Burlaga, 1974). We assme that the density and 
magnetic field intensity are constant at the inner boundary (O.l AU), 
being 800/cm^ and 5007 respectively. The results which follow refer to 
streams defined in this way, which we shall call our "standard" values. 

In general, both the direction and magnitude of ^ will depend on ^ 
at r = r^. Let us consider three cases: 1) the base of the field line 

in each volume element corotates with the sun, so that $Q=tan ^(B^/B^)=tan 
I 2) the field direction is the same for all volume elements at 
Tq, equal to 93°; 3) the field direction is constant, equal to 99^^ • 
such streams, the maximum value of the parameter fi ~ l/[l - (r-rj^)/RQ] 

= 3 (where r = 1 AU, r^ = 0.1 AU, = 1 AU/rad, and - kd’j/Y^ = ,75). 
The time profiles of V, n, B, and $ which an observer at 1 AU should see, 
according to our model, are shown in Figure 3. The speed profile is 
asymmetric because the fast elements overtake the slow elements ahead 
of them; the density is enhanced where V is increasing, with a maximum 
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value three times the ambient value; and the density is reduced in the 
region where V is decreasing. This is in quantitative agreement with 
observations and with the results of previously published, gas-dynamic 
models. In addition, our model predicts that the magnetic field inten- 
sity is enhanced in the region of increasing V, with a maximum a few 
times the ambient value at 1 AU; the magnetic field intensity is 
depressed below the ambient value at 1 AU in the region of decreasing V; 
and the magnetic field direction shows deviations from the spiral direc- 
tion which vary systematically with the speed profile. Similar results 
were obtained from the non-linear MHD models of Matsuda and Sakurai (1972) 
and Nakagawa and Wellck (1973) for a field which is parallel to the 
velocity in the corotating system. Our model shows quite clearly that 
these effects are kinematic rather than dynamical, and it has the 
advantage of allowing us to compute rather simply the effects of 
arbitrary field orientations near the inner boundary. 

Figure 3 shows that the magnitude and direction of B at 1 AU 


depends sensitively on at Tq- For Parker’s boundary condition ($^ = 
tan”^ ) one finds an enhancement of B in the region of increas- 

ing V, the maximum being«=4 times the value which would have occurred if 
a/v were zero (see Figure 3)* The value of corresponding to the mean 
speed (400 km/s) in this case is = ^ 6 ^ and BqCos 96 ° = BpSin 

(90° - Ij^) =« - 6*^ X (n/180) X Bq. If is constant equal to 93° , the 


enhancement is smaller because B^ is smaller by a factor of 2 (B^ 

- 3 ° X (ti/i8o) X Bq) while B « BQSin 96 ® « B does not change significantly. 


Similarly, if is constant equal to 99 ° » is larger because B^ i 
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larger while B^, remains approximately Bq. There is a corresponding 

change in angles » For Parker's boundary condition, one sees the angle 

change rapidly in the interaction region from a "tight spiral" where V 

is low to a "loose spiral" where V is high. For = 93'' a- similar 

change occurs, but in this case the field is generally more radial 

because of the smaller . Similarly, for = 95^ > 'the field is less 

radial because is larger, 

^o 

Observations show that the field is often at a large angle to the 
spiral direction. One possible interpretation is that it is the image of 
large deviations from the radial direction at the source. This leads us 
to study kinematic effects on magnetic fields of arbitrary initial 
orientation. We have calculated $ and jB|/n at 1 AU for a volume 
element starting at 0.1 AU as a function of the initial angle for 
different values of a/V. The results are shown in Figure 4 and 5 * 

Figure 4 shows $ at 1 AU as a function of at 0.1 AU. The 

dashed curve is the case a = 0, The variation of $ is due to the 

rotation of B with since Bj,^ = BoSin and B^ = B^cos (see 

(22)), When Bq = Bq r = 90° ) , ^n^ = ^ q /% ^ = 9CP ) » When 

= 96°, $Q = 138°, corresponding to the classical spiral result. When 

Bq = Bq ^ (^Q = 180°), 5 is also 180°, i.e,, an initially azimuthal field 

remains azimuthal if a/V = 0. As $q increases from 96° to l8o° , the 

field becomes more azimuthal because is increasing while B is 

Wo o 

decreasing. When a ^ 0, the shape of § ($q) is similar to that for a = 0, 
but it is displaced vertically and distorted somewhat (see Figure 4 ). 

This is most simply understood by considering ~ (^rQ “ * 


-13- 



Eqs. (22) and (24) show that when a/v > 0 the velocity shear twists the 
field so that there is a component |A|r which yields § < l8o° . When 
a/V < 0 the shear is in the opposite direction (-r) and $ > l8cf . 

The dependence of l^l/n = B/n on is shown for our standard 
stream in Figs. 5a and 5h. Note that B/n (i^ + l8cP ) = B/n ($q), 
meaning that B/n does not depend on the sense of ^ and we need only 
consider cP ^ s iSo*^ . When a = 0 (A = 0) , ^n^ = (B^/n^) sin r + 
(Bq/uq) ( r/r^) cos 0 and the variation of B/n is simply "sinusoidal" 

with a maximum at = OP ^ l8.0° . When A / 4 the position and magnitude 

of the maxima and minima change somewhat, as can be seen by considering 
(22). If one interprets visual coronal features as delineating magnetic 
field lines, it is to be expected that the field will be nearly radial 
most of the time so that values of close to 9CP and 270° in Figures 
4 and 5 are the most interesting. 



h, VAEIATIORS WITH DISTANCE FROM THE SUN 
a) Time Profiles at 0.3 AU 

Consider a stream which at 0.1 AU is characterized by V = + AV 

cos (46), V = 400 km/s, AV = 75 km/s, B = 5007, and n == 800 cm“^. 
o o o 

The time profiles of such a stream at one AU were shown in Figure 3 and 
were discussed in the preceding section. The time profiles of this 
stream at 0.3 AU are shown in Figure 6« The asymmetry of the speed 
profile is scarcely noticeable, but it is sufficient to produce a 30^ 
variation in the density. There is a mod\alation in B which depends on 
being <=« 15 ^ for $Q=tan"^ (-r^/R^), somewhat larger for a more azimu- 
thal initial field = 99^)? somewhat smaller for a more radial 

field (Iq = 93°). The variations in the direction of B at 0,3 AU also depend 
on When = tan”'^ (-r^/R^) the angle is more azimuthal at low 

speeds and more radial at high speeds. However, this is a small change 
and occurs across the interaction region, which is broad at 0.3 AU 
because the stream has not steepened much as yet. Consequently, the 
change in $ due to the stream is difficult to observe at 0.3 AU. The 
general conclusion from Figures 3 and 6 is that although there are 
qualitative similarities between the stream profiles at 0.3 and 1,0 AU, 
the quantitative differences are appreciable — the perturbations in n, B, 
and $ are relatively small at 0.3 AU but are large, obviously non-linear 
features at 1.0 AU. If large perturbations in B are observed near 
0.3 AU, it is safe to assume that they are not due to corotating streams, 
and alternative causes, such as variations in the source field or fields 
in flare-associated streams (Barouch et al., 1973), should be considered. 
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b) Variation of <By>' and mth r 


Radial variations of and in the solar wind have been reported 
by several observers. The custom is to compute and <B 0 > over =• 27 

day intervals for several different r and to compare the results with 

_p 

Parker’s model for a homogeneous solar wind which predicts Bj, ~ r and 
B 0 ~ r"^. The averaging approach and comparison with Parker’s model can 
be misleading, however, because of the presence of streams and the pos- 
sible variability of the field near the sun. 

Consider a solar wind made up of oar "standard" streams and consider 
sin 0=1. We have computed time averages of B 0 and B^ over these 
streams at several different distances, r, for = 93 °, 99 °, and Iq = 
tan”^ (-r^/R_). The results are shown in Figure 7. 

We find that <Bj.> r”^ with no noticeable dependence on 

Equation (22) shows that, this is because when $q is near 90 ^ Bj. *=» n^X 
Bro/rio = (’^oA)^^o ® ^ (^oA)^j where a = - 90° is a small 

angle. An r“^ dependence of <Bj,> has been reported by Burlaga and 
Hess ( 1968 ) and Villante and Mariani (1975) between 0.8 AU and 1 AU, by 
Behannon et al. (197^) between 0.^6 AU and 1 AU, and by Smith (1974) 
between 1 AU and 5 AU. (The Mariner 4 observations between 1 AU and 
1.5 AU (Coleman et al. , I 965 ) and the Mariner 5 observations between 
0.7 AU and 1 AU (Rosenberg, I 970 ; Rosenberg and Coleman, 1973) hid not 
give an r~^ dependence, however; the reason for this is not known). 

The behavior of <B 0 > is more ccmplicated than that of Our 

model shows that if $q is constant throughout the stream or equal to 
tan"^ (-r^/Rg), then <B 0 > ~ r"^, but it also shows that *=B 0 > is very 
sensitive to The reason is apparent from (22) which shows that 
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^0 " M Bq sin a (r^/r) ^ ^ B„ (vjr) a; i.e., B0 

varies as r“^, but it is also directly proportional to $q- 9CP . Thus, 
variations in are reflected proportionally in When §q takes 

values between 93^ and 95^ ? for oixr standard stream, <B0> can lie 
somewhere in the shaded area of Fig. 7* Measured values of <B0> depend 
on the value of near the sun, and on its fluctuations, as well as on 
the stream parameters. As the initial value and its statistical pro- 
perties may depend on time and on position, measurements of <B0>(r) per- 
formed during an extended period may well deviate significantly from an 

r dependence. Several groups have reported on the radial dependence 

-T 

of <B0^, expressed in the form <B0> <=: r . Villante and Mariani find 
7 = 2.5 + 2, Behannon et al. (197^^) find 7 = 1.22 or l.U, Smith (1974) 
finds results consistent with 7-4, and Rosenberg and Coleman found 
7 =1.85 aboard Mariner 5 and 1.22 aboard Mariner 4. These apparently 
conflicting results may well find their explanation in different distri- 
butions of during the time of these measurements. 

d) Map of Magnetic Field Intensity 
Figure 8 shows "equi- intensity" contours of B/B^ between 0.1 AU 
and 1.0 AU, computed for our standard stream (V = 400 + 75 cos 48) with 
§Q given by Parker's boundary condition. is the intensity that 

would be obtained if a = 0. The intensification of B in the region of 
increasing V is evident, as is the depression in the region of decreas- 
ing V. The depression is relatively small everywhere. The enhancement 
is large and increases rapidly near 1 AU. 
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Figure 8 was obtained by computing B (t) profiles at 20 different 

values of r, plotting ^ along circles with the various radii (assuming 

•°o 

corotation) , and connecting points with equal values of the ratio of 
Bfr)/Bo(r). There is some uncertainty in the shapes of the contours 
where the curvature is large, due to the finite mesh size that was 
used (0.05 AU), but generally, Figure 8 conveys a reasonably accurate 
image of the magnetic field intensity pattern in a corotating stream. 
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5. SUMVIARY MD DISCUSSION 


We have presented a theory for the kinematic hehavior of magnetic 
fields in streams in the solar wind between 0.1 AU and 1 AU, The theory 
is based on the assumption that V is constant and radial for any given 
volume element in the region considered and that the local speed gradients 
are not too large. This breaks down beyond 1 AU and for parts of the 
steepest streams at 1 AU, but generally it is a very good approximation 
for examining the behavior of ^ in the region inside 1.0 AU, By start- 
ing with speed profiles rather than temperature profiles, we fail to 
obtain well-known features concerning the temperature, but these have 
been adequately studied elsewhere and are not important as regards the 
kinematics of the magnetic field. 

The magnetic field variations are in the same sense as the density 
variations, but the ratio B/n depends strongly on the initial orientation, 
and the velocity gradient. The directional changes incurred by B are quite 
complex to describe, but are very sensitive to for near 90° • They 
are the result of the compression and "shear" of the component of B, 
which are caused by the velocity gradients of the stream. 

The implications of these calculations to the theory of cosmic ray 
propagations may be mentioned. Streams are a permanent feature of the 
interplanetary medium, and we expect regions of relatively intense magnetic 
field to be associated with each stream. These must be taken into account 
in cosmic ray propagation processes in a different manner than the usual 
diffusion -type calculations, as has already been pointed out, (Barouch 
and Raguideau, 1970; Barouch and Burlaga, 197^)* 
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Motivated by the availability of data obtained between 0,k6 and 
1 AU by the Mariner -Venus -Mercury spacecraft and the forthcoming data to 
0.3 AU from the Helios spacecraft, we have examined how the magnetic 
field should change in the ecliptic plane between 0.3 AU and 1 AU as a 
result of stream- kinematics. At 0.3 AU, the time variations are nearly 
linear, with small enhancements and depressions in B and n and with 
very small changes in $ across the stream. These variations grow non- 
linearly as one approaches 1 AU. The growth depends on the field 
direction, at the inner boundary. It has been the practice to 
compute time averages of and B 0 at different distances and compare 
them with the r and r"*^ variations predicted by the classical spiral 
model for a homogeneous wind with no streams. Averaging over a model 

_P 

stream, we find that <Bj.> ~ r and is insensitive to $q, in agreement 
with the spiral model and with some observations; but we find that 
<B 0 > is very sensitive to and will not generally be proportional to 
r"^ if changes in the stream, which may explain why the observations 
do not show an r“^ dependence and vary among themselves. The field direc 
tion is very sensitive to the orientation near the sun and small depar- 
tures from the "spiral angle", $o = tan"^ (-ro/R^), cause a large spread 
about the spiral angle at r ^ r^, the spread increasing with r . Thus, 
the observed variation in the direction of B at 1 AU might be partly 
the result of relatively small fluctuations in the direction of B near 
the sun . 
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FIGURES 


1, This illustrates the geometry and defines angles. Note that V(8) is 
stationary in the corotating frame. 

2, (top) fx. (a/v) at 0.3 AU, 0,5 AU, and 1.5 AU (bottom) (X (R) between 

0.1 AU and 1 AU for a/v = 0, + 0,4, + 0.7- /u. is n(r)/n^ where n^ 

is the density at 0.1 AU. 

3, Time profile of a stream at 1 AU. The magnetic field intensity, B, 
and direction $, depend on the angle at r = 0.1 AU, but n is 
independent of 

4, $ at 1 AU(ij_) versus at 0.1 AU for a/v = 0, + 0,4, + 0.7. 

5, B/n versus a) a/V > 0, b) a/v < 0. 

6, Time profile of a stream at 0.3 AU. 

7, an averaged over the time profile of a "standard" stream 

as a function of distance from the sun, E, ~ R“^j but <B0(R)> 

depends on 0^. 

8, Contour map of B/Bg^_Q for our "standard" stream. Bg^^ is the value 
of B(r,0) that would be measured in the absence of a stream. The 
view shown is in the ecliptic plane. 
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APPENDIX 


Validity of the Kinematic Approximation 


All of the results discussed above are based on the kinematic 
approximation, in which we set the right hand side of the radial com- 
ponent of the momentum equation 

J olt 1)'^ 

equal to zero. Neglecting the (jg,*V) term is justifiable because 

it is on the order of the ratio of the size of the interaction fegion to 

the radius of curvature of the spiral field line, which is < 0.1. Since 

we are interested only in an order of magnitude estimate of the ratio of 

the RHS to the LHS of (A.l), we can make the following approximations to 

estimate the size of the RHS of (A.l); B = 

2 

and P = Pj^ (where P is the pressure) which is an adiabatic 

law with 7=2. Let us divide both sides of A.l by ^ and set ~ V — , 

dt ^ 

which is valid when V = constant* The RHS 

-V 2 

is then — JL where V = V ^ + V. ^ = (2P, /Q ) + (B.^/(h-nQ ’I 'i at 


is then 


^ g, where V = t = (2P^/e^) t (B,_ /(Un?^)) at 

X XIX 


0*1 AU; here, and V are the sound speed and AlfvSn speed respectively. 

S A 

2 2 2 

The equation can then be integrated to give - V -V 2(1 - ), 

O O 1 

where the subscript zero refers to the inner boundary. Since = 

1 ^1 
is the magnetoacoustic speed at 1 AU in the absence of distortions due 

to streams (a = o) . Then 


(vo% a -3.^ Vrti 
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The second terra on the RHS of (A. 2) is the result of an acceleration of 

the volurae element by the ambient gradients of p and B, and the third 

term represents the effect of a stream. ¥e consider that the effect of 

2 2 

a stream is negligible if % ^1 < 0 . 0<V^ . For the parameters that we 

•2 2 2 S 2 

have been using, ^ 2 x 10^ (km/s) and «« 2 x 10^, (km/s) , so 

the effects of streams on V (the radial component of the velocity) can 

indeed be neglected to first approximation up to 1 AU if | ^| < 0.7. This 

is true for most streams < 1.0 AU, but the approximation probably breaks 

down for the steepest streams at 1 AU, and it cannot be used much beyond 

1 AU for most streams. 

It is interesting to note from (A. 2) that when our approximation 

does break down, energy is extracted from the flow to provide the potential 

energy in the pressure pulse, the decrease in V being greatest where 

i.e. the density enhancement, is greatest. This implies that the speed 

profile which is computed assuming ~ ~ 0 will be altered such that it 

dt 

is steeper at 1 AU, and the density enhancement will appear to be closer 


to the front of the resulting stream. 
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